Transition to Instability in a Periodically Kicked Bose-Einstein Condensate on a Ring 
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A periodically kicked ring of a Bose-Einstein condensate is considered as a nonlinear generalization 
of the quantum kicked rotor, where the nonlinearity stems from the mean field interactions between 
the condensed atoms. For weak interactions, periodic motion (anti-resonance) becomes quasiperiodic 
(quantum beating) but remains stable. There exists a critical strength of interactions beyond which 
quasiperiodic motion becomes chaotic, resulting in an instability of the condensate manifested by 
exponential growth in the number of noncondensed atoms. In the stable regime, the system remains 
predominantly in the two lowest energy states and may be mapped onto a spin model, from which 
we obtain an analytic expression for the beat frequency and discuss the route to instability. We 
numerically explore parameter regime for the occurrence of instability and reveal the characteristic 
density profile for both condensed and noncondensed atoms. The Arnold diffusion to higher energy 
levels is found to be responsible for the transition to instability. Similar behavior is observed for 
dynamically localized states (essentially quasiperiodic motions), where stability remains for weak 
interactions but is destroyed by strong interactions. 

PACS numbers: 03.75.-b, 05.45.-a, 03.65.Bz, 42.50.Vk 
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I. INTRODUCTION 



The (5-kicked rotor is a textbook paradigm for the 
study of classical and quantum chaos Qj. In the classi- 
cal regime, increasing kick strengths destroy regular pe- 
riodic or quasiperiodic motions of the rotor and lead to 
the transition to chaotic motions, characterizing by diffu- 
sive growth in the kinetic energy. In quantum mechanics, 
chaos is no longer possible because of the linearity of the 
Schrodinger equation and the motion becomes periodic 
(anti- resonance) , quasiperiodic (dynamical localization), 
or resonant (quantum resonance) 0, Q • Experimental 
study of these quantum phenomena have been done with 
ultra-cold atoms in periodically pulsed optical lattices 
0. However, most of previous investigations have been 
focused on single particle systems and the effects of in- 
teraction between particles have not received much at- 
tention dEl. 

In recently years, the realization of Bose-Einstein con- 
densation (BEC) 0| of dilute gases has opened new op- 
portunities for studying dynamical systems in the pres- 
ence of many-body interactions. One can not only pre- 
pare initial states with unprecedented precision and pure- 
ness, but also has the freedom of introducing interactions 
between the particles in a controlled manner. A natural 
question to ask is how the physics of the quantum kicked 
rotor is modified by the interactions. In the mean field 
approximation, many-body interactions in BEC are rep- 
resented by adding a nonlinear term in the Schrodinger 
equation y| (such nonlinear Schrodinger equation also 
appears in the context of nonlinear optics |9j ) . This non- 
linearity makes it possible to bring chaos back into the 
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system, leading to instability (in the sense of exponential 
sensitivit y to initial conditions) of the condensate wave 
function • The onset of instability of the condensate 
can cause rapid proliferation of thermal particles that 
can be observed in experiments 12]. It is therefore im- 
portant to understand the route to chaos with increasing 
interactions. This problem has recently been studied for 
the kicked BEC in a harmonic oscillator 

In Ref. |l3l |. we have investigated the quantum dy- 
namics of a BEC with repulsive interaction that is con- 
fined on a ring and kicked periodically. This system is 
a nonlinear generalization of the quantum kicked rotor. 
From the point of view of dynamical theory, the kicked 
rotor is more generic than the kicked harmonic oscilla- 
tor, because it is a typical low dimensional system that 
obeys the KAM theorem, while the kicked harmonic os- 
cillator is known to be a special degenerate system out 
of the framework of the KAM theorem ^ij- It is very 
interesting to understand how both quantum mechanics 
and mean field interaction affect the dynamics of such a 
generic system. 

In this paper, we extend the results of Ref. in- 
cluding a more detailed analysis of the model consid- 
ered there as well as new phenonmena. We will focus 
our attention on the relatively simpler case of quantum 
anti-resonance, and show how the state is driven towards 
chaos or instability by the mean field interaction. The pa- 
per is orgnaized as follows: Section II lays out our physi- 
cal model and its experimental realization. Section III is 
devoted to the case of weak interactions between atoms. 
We find that weak interactions make the periodic motion 
(anti-resonance) quasiperiodic in the form of quantum 
beating. However, the system remains predominantly in 
the lowest two energy levels and analytic expressions for 
the beating frequencies are obtained by mapping the sys- 
tem onto a spin model. Through varying the kick period, 
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we find the phenomenon of anti-resonance may be recov- 
ered even in the presence of interactions. The decoher- 
ence effects due to thermal noise are discussed. Section 
IV is devoted to the case of strong interactions. It is 
found that there exists a critical strength of interactions 
beyond which quasiperiodic motion (quantum beating) 
are destroyed, resulting in a transition to instability of 
the condensate characterized by an exponential growth 
in the number of noncondensed atoms. Universal criti- 
cal behavior for the transition is found. We show that 
the occurrence of instability corresponds to the process of 
Arnold diffusion, through which the state can penetrate 
through the KAM tori and escape to high energy levels 
[Tfij . We study nonlinear effects on dynamically localized 
states that may be regarded as quasiperiodic . Similar 
results are obtained in that localization remians for suf- 
ficiently weak interactions but become unstable beyond 
a critical strength of interactions. Section V consists of 
conclusions. 



billiards 01 > or optical-dipole traps produced by red- 
detuned Laguerre-Gaussian laser beams of varying az- 
imuthal mode index |l9j . The <5-kick may be realized 
by adding potential points along the ring with an off- 
resonant laser £|, or by illuminating the BEC with a 
periodically pulsed strongly detuned running wave laser 
in the lateral direction whose intensity / is engineered 
to I = Iq + Iix/R 20]. In the experiment, we have the 
freedom to replace the periodic modulation cos (9) of the 
kick potential with cos(j0), where j is an integer. Such 
replacement revises the scaled interaction constant g and 
kicked strength K, but will not affect the results obtained 
in the paper. The interaction strength g may be adjusted 
using a magnetic field-dependent Feshbach resonance or 
the variation of the number of atoms |2l| ■ 



B. Dynamics of non-condensed atoms: Bogoliubov 
equation 



II. PHYSICAL MODEL: KICKED BEC ON A 
RING 

Although the results obtained in this paper are com- 
mon properties of systems whose dynamics are governed 
by the nonlinear Schrodinger equation, we choose to 
present them here for a concrete physical model: a kicked 
BEC confined on a ring trap, where the physical mean- 
ings of the results are easy to understand. The descrip- 
tion of the dynamics of this system may be divided into 
two parts: condensed atoms and non-condensed atoms. 

A. Dynamics of condensed atoms: Gross-Pitaevskii 
equation 

Consider condensed atoms confined in a toroidal trap 
of radius R and thickness r, where r <C R so that lat- 
eral motion is negligible and the system is essentially 
one-dimensional [17|. The dynamics of the BEC is de- 
scribed by the dimensionless nonlinear Gross-Pitaveskii 
(GP) equation, 

i§^= (-\^+9\i>? + K™s68 t {T)\i>, (1) 

where g = 8NaR/r 2 is the scaled strength of nonlinear 
interaction, N is the number of atoms, a is the s-wave 
scattering length, K is the kick strength, 8t(T) represents 
^2 6 (t — nT), T is the kick period, and 9 denotes the az- 

n 

imuthal angle. The length and the energy are measured 
in units R and -^s- 5 respectively. The wavefunction 

ij)(6,t) has the normalization f ( 2ir \if;\ 2 d9 = 1 and satis- 
fies periodic boundary condition ip (9, t) = ip (9 + 2ir, t). 

Experimentally, the ring-shape potential may be real- 
ized using two 2D circular "optical billiards" with the 
lateral dimension being confined by two plane optical 



The deviation from the condensate wave function is de- 
scribed by Bogoliubov equation that is obtained from a 
linearization around GP equation [Til ]. In Castin and 
Dum's formalism, the mean number of noncondensed 
atoms at zero temperature is described by (5N(t)) — 
Y,kLi( v k(t)\ v k(t)) , where Ufc(i) are governed by 

4{::)={-% -%){::)■ < 2 » 

where Hi = p 2 /2 + g\iP\ 2 - (i + gQ\ip\ 2 Q + Kcos65 t (T), 
H-2 = gQip 2 Q*, M is the chemical potential of the ground 
state, i/j is the ground state of GP equation and the pro- 
jection operators Q are given by Q = 1 — \ip)(tp\- 

The number of noncondensed atoms describes the de- 
viation from condensate wavefunction and its growth rate 
characterizes the stability of the condensate. If the mo- 
tion of the condensate is stable, the number of noncon- 
densed atoms grows at most polynomially with time and 
no fast depletion of the condensate is expected. In con- 
trast, if the motion of the condensate is chaotic, the num- 
ber of noncendensed atoms diverges exponentially with 
time and the condesate may be destroyed in a short time. 
Therefore the rate of growth of the noncondensed atoms 
number is similar the Lyapunov exponent for the diver- 
gence of trajectories in phase space for classical systems 



III. WEAK INTERACTIONS: 
ANTI-RESONANCE AND QUANTUM BEATING 

In this section, we focus our attention on the case of 
quantum anti-resonance, and show how weak interactions 
between atoms modify the dynamics of the condensate. 
Quantum anti-resonace is a single particle phenomenon 
characterized by periodic recurrence between two differ- 
ent states, and its dynamics may be described by Eq. (JTJ 
with parameters g = and T = 2n Q . 
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Figure 1: Plots of average energy E(t) versus the number of 
kicks t for three values of g. The kick strength K — 0.8. 



Figure 2: Fourier transformation of the energy evolutions in 
Fig.l. The Unit for the frequency is 1/T. 



A. Quantum Beating 

In a non-interacting gas, the energy of each particle 
oscillates between two values because of the periodic re- 
currence of the quantum states. In the presence of inter- 
actions, single particle energy loses its meaning and we 
may evaluate the mean energy of each particle 



25(f) = 



2 IT 



(3) 



To determine the evolution of the energy, we numeri- 
cally integrate Eq.Q over a time span of 100 kicks, using 
a split-operator method |22j . with the initial wavefunc- 
tion ip being the ground state ip(9, 0) = l/^/2n. After 
each kick, the energy 25(f) is calculated and plotted as 
shown in Fig.l. 

In the case of non-interaction (Fig. 1(a) g = 0), we 
see that the energy 25(f) oscillates between two values 
and the oscillation period is 2T, indicating the periodic 
recurrence between two states (anti- resonance) . 

The energy oscillation with weak interaction (g = 0.1) 
in Fig. 1(b) shows a remarkable difference from that for 
non-interaction case. We see that the amplitude of the 
oscillation decreases gradually to zero and then revives, 
similar to the phenomena of beating in classical waves. 
Clearly, it is the interactions between atoms in BEC that 
modulate the energy oscillation and produce the phenom- 
ena of quantum beating. As we know from classical waves, 
there must be two frequencies, oscillation and beat, to 
create a beating. These two frequencies are clearly seen 
in Fig. 2 that is obtained through Fourier Transform of 
the energy evolutions in Fig.l. For the non-interaction 
case (Fig. 2(a)), there is only the oscillation frequency 



lose — 0.5 /T, corresponding to one oscillation in two 
kicks. The interactions between atoms develop a new 
beat frequency fbeau a $ well as modify the oscillation 
frequency f 08C , as shown in Fig. 2(b). 

In Fig. 3, these two frequencies are plotted with re- 
spect to different interacting constant g and kick strength 
K. We see that both beat and oscillation frequencies vary 
near linearly with reapect to the interaction strength g. 
More interestingly, the two frequencies satisfy a conser- 
vation relation 



fosc + fbeat/% — 1/2- 



(4) 



For a strong interaction [Fig. 1(c)], i.e., g > 1.96, we 
find that the energy's evolution demonstrates an irregular 
pattern, clearly indicating the collapse of the quasiperi- 
odic motion and the occurrence of instability. The corre- 
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Figure 3: Plots of beat and oscillation frequencies versus the 
interaction strength (a) and kick strength (b), where the scat- 
ters are the results from numerical simulation using GP equa- 
tion and lines from analytic expression Eqs. il l ■">> and llfi^ . 



4 



sponding Fourier transformation of the energy evolution 
(Fig. 2(c)) has no sharp peak. This transition to instabil- 
ity will be discussed in Section IV in details. 



B. Spin model 

The phenomena of quantum beating can be understood 
by considering a two-mode approximation [§3 to the GP 
equation. In this approximation, condensed atoms can 
only effectively populate the two lowest second quan- 
tized energy modes. The validity of this two-mode model 
is justified by the following facts which are observed in 
the numerical simulation. First, the total energy of the 
condensate is quite small so that we can neglect the pop- 
ulation at high energy modes and keep only the states 
with quantized momentums and ±1. Second, the to- 
tal momentum of the condensate is conserved during the 
evolution. Therefore the populations of the states with 
momentum ±1 are same if the initial momentum of the 
condensate is zero (the ground state). 

Here we consider a quantum approach of this two-mode 
model, which yields an effective spin Hamiltonian in the 
mean field approximation. By considering the conserva- 
tion of parity we may write the Boson creation operator 
for the condensate as 



1> 



t _ 



/2?r V 



(5) 



where cv and w are the creation operators for the ground 
state and the first excited states, satisfying the commuta- 
tion relation [fi, fiT] = 1, b,b' = 1, and particle number 

conservation a^a + wb — N . 

Substituting Eq. JSJ into the many-body Hamiltonian 
of the system 



H = 



2tt 



+V> f (K cos 9S T (t))^ 
we obtain a quantized two-mode Hamiltonian 



(6) 



H ~2 Lz ' AttN 



-V2KLJ t (f) 



(7) 



in terms of the Bloch representation by defining the an- 
gular momentum operators, 



L, 



tfb + ati 



a+S — aSt 



(8) 



a^a-tib 



where we have discarded all c-number terms. 

The Heisenberg equations of motion for the three an- 
gular momentum operators reads 



dL x 




~dT 


2 


dL y 


L x 


~df 


2~ 



87riV 
V2KL Z S T (t) 



-L v - JL- (jL, L z ] (9) 

47V v 8ttN X v J V ; 

{/.../.,} 

87TiV I ' I 



-L x 



= ^{L y ,L x } + V2KL v 6 T (t) 

where = LiLj + LjLi. 

The mean field equations for the first order expec- 
tion values of the angular momentum operators 

are obtained by approximating second order expectation 
values (^LiLj^j as products of (i>i) an d (^}' that 

(LiLj^j = (j^i) (^j) Defining the single-particle 

Bloch vector 



S = (S x , S y , S z ) = - \^L X ) , (L v ) , (L z ) j , (10) 
we obtain the nonlinear Bloch equations 

S. = |± + £-£s,U, (11) 



SyS x + V2KS y S T (t) , 



where we have used N 3> 1. The mean field Hamiltonian 
in the spin representation reads 



s, si 

4 + 8 



V2KS x S t (T). (12) 



From the definition of the Bloch vector, we see 
that S z corresponds to the population difference and 
— axctan(S v / S x ) corresponds to the relative phase a be- 
tween the two modes. This Hamiltonian is similar to a 
kicked top model 0, but here the evolution between 
two kicks is more complicated. 

With the spin model, we can readily study the dynam- 
ics of the system. For the non-interaction case (g = 0), 
The Bloch equations (|11|) become 



S x — ^S y , 



(13) 



S y = ~S x -V2KS z 6 T (t), 

S z = V2KS y 6 T (t). 

We see that the evolution between two consecutive kicks 
is simply an angle tt rotation about the z axis, which 
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Figure 4: Periodic stroboscopic plots of the projection of spin 
on the S z — plane. The thick line and dots correspond to 
the orbits with initial spin S = (0, 0, 1). K = 0.8 (a) g = 0.0; 
(b) ff = 0.1. 

yields the spin transformation S x — > — S x , S y — > — S y . 
The spin initially directing to north pole (S = (0,0, 1)) 
stays there for time duration T, then the first kick rotates 
the spin by an angle y/2K about the x axis and now S = 
(0, - sin (V2K) , cos (V2K)). The following free evolu- 
tion rotates the spin to S = (0, sin [\/2K^ , cos (\/2-K^) ■ 
Then, the second kick will drive the spin back to north 
pole through another rotation of y/2K about the x axis. 
With this the spin's motion is two kick period recurrence 
and quantum anti-resonance occurs. 

The motion of the spin is more complicate with interac- 
tions. The spin co mponen ts at x and y directions can be 
written as S x = y/l — S% cos a and S y = — ^1 — S z sin a 
in terms of population difference S z and relative phase 
a, which yields the relation 

• g 2 

S x = aS y — —S Z S V cos a 

TT 

during the free evolution. Comparing this equation with 
the Bloch equation , we obtain the equation of mo- 
tion for the relative phase 

We see the motion between two consecutive kicks is ap- 
proximately described by a rotation of 7t + g(l + 3S z )/4 
about the z axis. Compared with the noninteraction case, 
the mean field interaction leads to an additional phase 
shift g(l + 3S z )/A. This phase shift results in a deviation 
of the spin from S x = plane at time 2T _ , i.e., moment 
just before the second kick. As a result, the second kick 
cannot drive the spin back to its initial position and quan- 
tum anti-resonance is absent. However, the phase shift 
will be accumulated in future evolution and the spin may 
reach S x = plane at a certain time mT~ (beat period) 
when the total accumulated phase shift is ir/2. Then the 
next kick will drive spin back north pole by applying an 
angle y/2K rotation about the x axis. 
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Figure 5: (a) Plots of relative phase versus the number of 
kicks t, where g = 0.1, K = 0.8. (b) Schematic plot of the 
phase shift. nT"' + ' represents the moment just before (after) 
the nth kick. 



The above picture is confirmed by our numerical so- 
lution of the spin Hamiltonian with fourth-order Runge- 
Kutta method [25l [. In Fig. 4, we plot the phase por- 
traits of the spin evolution by choosing different initial 
conditions of the population difference S z and relative 
phase a. Just after each kick three spin components 
Si (i = x, y, z) are determined and their projections on 
S z = plane are plotted. For the noninteraction case 
(Fig. 4(a)), we see the motion of the spin is an osillation 
between the north pole A and another point -B, indi- 
cating the occurence of quantum anti-resonance. The 
interaction between atoms changes the phase portraits 
drametically (Fig. 4(b)). Around the north pole, a fixed 
point surrounded by periodic elliptic orbits appears. The 
two-point oscillation is shifted slowly and form a contin- 
uous and closed orbit, representing the phenomenon of 
quantum beating. 

In Fig. 5 we see that the relative phase at the moment 
just before the even kicks increases almost linearly and 
reaches 2ir in a beat period. The slope of the increment 
reads, j RP = (a (4T~) - a (2T~)) /2, which can be de- 
duced analytically. With this and through a lengthy de- 
duction, we obtain an analytic expression for the beat 
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Figure 6: (a) Plot of population difference versus the number 
of kicks t, where g — 0.1, K = 0.8. (b) Details of (a) at the 
middle of a beat period. The population difference decreases 
in two consecutive kicks. 
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frequency to first order in g, 



fbe 



9_ 

4tt 



(l + 3 cos 



(15) 



In Fig. 6, we plot the evolution of population differ- 
ence and the phenomenon of quantum beating is clearly 
seen in the spin model. Notice that there is one peak miss 
of the oscillation at the middle of one period because the 
population difference decreases in two consecutive kicks. 
Taking account of this missed peak, we obtain the oscil- 
lation frequency 



fo 



N, 



total 



N„ 



2N, 



total 



-i-if 
2 2 J fr ea * ' 



(16) 



where N to tai and N, 



are total and missed numbers of 
peaks, respectively. 

The analytical expressions Eqs. (|15|) and l|16[> of the 
beat and oscillation frequencies are in very good agree- 
ment with the numerical results obtained from G-P equa- 
tion, as shown in Fig. 3. Therefore the beating provides 
a method to measure interaction strength in an experi- 
ment. 
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Figure 8: Plots of population difference with respect to the 
number of kicks in the presence of decoherence. K — 0.8, g = 
0.2. Thin lines correspond to F — 0. Thick lines correspond 
to (a) 2TrF/g = 0.001, (b) 27vT/g = 1. 



In the spin model, we see that it is the additional phase 
shift originating from weak interactions that destroys 
the two kick period recurrence of the anti-resonance and 
leads to the phenomenon of quantum beating. There- 
fore we will still be able to observe the quantum anti- 
resonance even in the presence of interactions if the ad- 
ditional phase shift may be compensated. Actually, the 
additional phase shift can be cancealed by varying the 
kick period T so that the relative phase a only change tt 
between two consecutive kicks. Using Eq. <|14[) . we nn d 
the new kick period for anti-resonance in the presence of 
interactions may be approximated as 



Tar 



8tt 2 



4-7T - 



3g cos (\2K 



(17) 
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Figure 7: Plots of average energy E (t) versus the number of 
kicks t for K = 0.1, g = 0.1. Tar is determined through Eq.( 



In Fig. 7, we plot the evolution of the average energy 
E (t) with the new kick period Tar- We see that the 
energy oscillates between two values and the oscillation 
period is 2Tar, clearly indicating the recovery of the anti- 



resonance. 



D. Decoherence due to thermal noise 

In realistic experiments, the decoherence effects always 
exist. Generally, decoherence originates in the coupling 
to a bath of unobserved degree of freedom, or the inter- 
particle entanglement process j^]. The main source of 
decoherence in a BEC is the thermal cloud of particles 
surrounding the condensate. Thermal particles scatter- 
ing off the condensate will cause phase diffusion at a rate 
proportional to the thermal cloud temperature. Here we 
consider a simple model that accounts for the effect of 
the thermal noise on the two-mode dynamics by adding 
a r = 1/r transversal relaxation term 26] into the mean- 
held equations of the motion 



St. — 



Sy — 



2 8tt 8tt ' v 



1 



JL -1l 

8tt 8tt 



(18) 

s)j s x - V2ks z s t (t) - rs yi 



S z — —S y S x + V^KSySr (t) ■ 

IT 

In Fig. 8, we plot the evolution of the population dif- 
ference S z for different decoherence constant T. We see 
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that the phenomenon of quantum beating is destroyed by 
strong thermal noise (Fig. 8(b)), while survives in weak 
noise (Fig. 8 (a)). For large decoherence constant, the 
population difference S z decays to exponentially and 
the characteristic time is the just the decoherence time 
t = 1/r. Therefore the decoherence time r must be 
much larger than the beat period 2ir/ fbeat to observe the 
phenomenon of quantum beating, which yields 



2nT/g«±-(l 



3cos(V2iT)). 



(19) 



In the case of K = 0.8, Eq. H19|) gives an estimation 
2ttT I g << 0.2, which agrees with the numerical results 
shown in Fig. 8. 



IV. STRONG INTERACTIONS: TRANSITION 
TO INSTABILITY 

A. Characterization of the instability: Bogoliubov 
excitation 

Tuning the interaction strength still larger means en- 
hancing further the nonlinearity of the system. From 
our general understanding of nonlinear systems, we ex- 
pect that the solution will be driven towards chaos, in 
the sense of exponential sensitivity to initial condition 
and random evolution in the temporal domain. The lat- 
ter character has been clearly displayed by the irregular 
pattern of the energy evolution in Fig. 1(c). On the other 
hand, the onset of instability (or chaotic motion) of the 
condensate is accompanied with the rapid proliferation 
of thermal particles. Within the formalism of Castin and 
Dum ^l| described in Section II, the growth of the num- 
ber of the noncondensed atom will be exponential, simi- 
lar to the exponential divergence of nearby trajectories in 
phase space of classical system. The growth rate of the 
noncondensed atoms is similar to the Lyapounov expo- 
nent, turning from zero to nonzero as instability occurs. 

We numerically integrate Bogoliubov equation J5J for 
the Wfc, Vk pairs over a time span of 100 kicks, using a 
split operator method, parallel to numerical integration 
of GP equation The initial conditions 



«fc(0) 
Vfe(0) 



2 U-r 1 



ikO 



(20) 



for initial ground state wavefunction ip(0) = l/v2~7r, are 
obtained by diagonalizing the linear operator in Eq.Q 

El, where C=( fc2/2 fc ;g^ ) 1/4 . 

After each kick the mean number of noncondensed 
atoms is calculated and plotted versus time in Fig. 9(a). 
We find that there exists a critical value for the interac- 
tion strength, i.e., g c = 1.96, above which, the mean num- 
ber of noncondensed atoms increases exponentially, indi- 
cating the instability of BEC. Below the critical point, 
the mean number of noncondensed atoms increases poly- 
nomially. As the nonlinear parameter crosses over the 
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Figure 9: (a) Semilog plot of the mean number of noncon- 
densed atoms versus the number of kicks t. The thicker lines 
are fitting functions. K = 0.8, g = 0.1 (dashed line, fit- 
ting function 0.0003i 13 ), g = 1.5 (dotted line, fitting func- 
tion 0.001 It 2 ), g — 2.0 (dash dotted line, fitting function 
0.32exp(0.1t)). The inset shows the interaction dependence of 
the growth rate. The scatters are from numerical simulation 
and the solid line is the fitting function 0.33(g — t.96) 1 ^ 2 . (b) 
Phase diagram of the transition to instability. 



critical point, the growth rate turns from zero to nonzero, 
following a square-root law (inset in Fig. 9(a)). This scal- 
ing law may be universal for Bogoliubov excitation as 
confirmed by recent experiments [l2| 

The critical value of the interaction strength depends 
on the kick strength. For very small kick strength, the 
critical interaction is expected to be large, because the 
ground state of the ring-shape BEC with repulsive inter- 
action is dynamically stable 29]. For large kick strength, 
to induce chaos, the interaction strength must be large 
enough to compete with the external kick potential. So, 
in the parameter plane of {g, k), the boundary of insta- 
bility shows a "U" type curve (Fig. 9(b)). 




Figure 10: Plots of condensate and noncondensate densities, 
where K = 0.8. (a,b) g = 0.1; (c,d) g = 2.0. 
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Across the critical point, the density profiles of both 
condensed and noncondensed atoms change dramatically. 
In Fig. 10, we plot the temporal evolution of the den- 
sity distributions of condensed atoms as well as noncon- 
densed atoms. In the stable regime, the condensate den- 
sity oscillates regularly with time and shows clear beat- 
ing pattern (Fig. 10(a)), whereas the density of the non- 
condensed atoms grows slowly and shows main peaks 
around 9 = ±7r and 0, besides some small oscillations 
(Fig. 10(b)). In the unstable regime, the temporal oscil- 
lation of the condensate density is irregular (Fig. 10(c)), 
whereas the density of noncondensed atoms grows explo- 
sively with the main concentration peaks at 9 = ±ir/2 
where the gradient density of the condensed part is max- 
imum (Fig. 10(d)). Moreover, our numerical explorations 
show that the cos 2 9 mode (FiglO.(b)) dominates the 
density distribution of the noncondensed atoms as the 
interaction strength is less than 1.8. Thereafter, the 
sin 2 9 mode grows while cos 2 9 mode decays, and finally 
sin 2 9 mode become dominating in the density distribu- 
tion of noncondensed atoms above the transition point 
(FiglO.(d)). Since the density distribution can be mea- 
sured in experiment, this effect can be used to identify 
the transition to instability. 



B. Arnold Diffusion 

We have seen that strong interactions destroy beat- 
ing solution of the GP equation and the motion of the 
condensate become chaotic, characterized by exponential 
growth in the number noncondensed atoms. The remain- 
ing question is how the motion of condensate is driven to 
chaos, that is, the route of the transition to instability. 

The transition to chaos for the motion of the conden- 
sate can be clearly seen in the periodic stroboscopic plots 
of the trajectories for both two-mode approximation 
(Fig.ll(a-d)) and four-mode approximation (Fig.ll(e-f), 
it is exact for the interaction region we consider in Fig. 11) 
to the original GP equation (JJJ. The solution oscillates 
between two points A and B (point C is identical to 
A in spin model.) for noninteraction case, and forms 
closed path in the phase space for the weak interaction 
(Figll.(a)). Note that compared with the two- mode cal- 
culation, the effective interaction strength in the four- 
mode approximation is rescaled, to give the comparable 
pattern in the phase space. 

With increasing interaction strength, the stable quasi- 
periodic orbits in Fig. 11(b) bifurcates into three closed 
loops (Fig. 11(c)) and chaos appears in the neighborhood 
of the hyperbolic fixed points. However, diffusion from 
one stochastic region to another are still blocked by KAM 
tori for the two-mode system. In Fig. 11(d), the trajec- 
tory with initial condition S z — 1, a — is closer to the 
chaotic region. 

The above discussion is based on two-mode approxima- 
tion, actually, the solution is coupled with other modes 
of higher energy states. For small interaction, this cou- 



pling is negligible and the four-mode simulation gives 
the same results as seen in Fig.ll(a,b,e,f). For large 
interaction, this coupling is important and our system 
is essentially high-dimensional (d > 2). One important 
character of a high-dimensional dynamical system is that 
KAM tori (d-dimension) can not seperate phase space 
(2d-dimcnsion) and the whole chaotic region is intercon- 
nected. If a trajectory lies in a chaotic region it can 
circumvent KAM tori and diffuse to higher energy states 
through Arnold diffusion. This process is clearly demon- 
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Figure 11: Periodic stroboscopic plots of population difference 
with respect to the relative phase between the first two modes. 
K = 0.8. (a-d) the two-mode model, where (a) g = 0.1, 
circle dots corresponds to g — 0; (b) g = 1.5; (c) g — 1.9; 
(d) g — 2.0; The thicker line and larger dots on the phase 
portaits represent the trajectories with initial conditions S z = 
1, a = 0. (e-h) the four-mode approximation, the portait is 
the projection on the first two modes of the trajectory with 
initial condition S z = 1, a = 0. (e) g — 0.1; (f) g = 1.6; (g) 
g = 2.2; (h) g = 2.5. 
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strated in Fig.ll(g,h). We see that the trajectory diffuses 
along the separatrix layers, circumvents the KAM tori, 
and finally spreads over whole phase space (Fig. 11(h), 
g = 2.5). We also calculate the diffusion coefficient 

r> - 2 \ Em Z g "| 2 fo-n 

Ve ~ "777 — 77 "77 T' V 21 ) 

J(J — 1) — ' 1 (m — n) 

where E m is the energy after the mth kick, J is the total 
number of kicks. For g — 2.2 and g — 2.5, the diffusion 
rates are 7.2 x 1CP 11 and 1.5 x 1CP 9 respectively. 

Arnold diffusion allows the state to diffuse into higher 
energy states, which destroys quasiperiodic motion of the 
quantum beating and leads to the transition to instabil- 
ity. As Arnold diffusion occurs, the motion of the con- 
densate becomes unstable and the number of the non- 
condensed atoms grow exponentially, as we have seen in 
above discussion. 

Arnold diffusion is a general property of the nonlinear 
Schrodinger equation in the presence of strong interac- 
tions. However, it may be hard to observe the whole 
process of Arnold diffusion in realistic BEC experiments 
because of the limit number of atoms (10 6 ). As Arnold 
diffusion occurs, the instability of the condensate leads to 
the exponential growth of thermal atoms which destroy 
the condensate, as well as invalidate the GP equation JTJ 
in a short time; while clear signature of the whole Arnold 
diffusion process may only be observed in a realtive long 
period. On the other hand, Arnold diffusion may be 
observed in the context of nonlinear optics, where the 
GP equation describes the propagation of photons. 
The number of photons is very large and the interactions 
between them are very weak, therefore it is possible to 
have a long diffusion process without invalidating the GP 
equation QJ. 

C. Dynamical localized states 

Although the above discussions have been focused on 
a periodic state of anti-resonance, the transition to in- 
stability due to strong interactions also follows a similar 
path for a dynamically localized state. The only differ- 
ence is that we start out with a quasiperiodic rather than 
periodic motion in the absence of interaction. This means 
that it will generally be easier to induce instability but 
still requires a finite strength of interaction. 

In Fig. 12, we show the nonlinear effect on a dynami- 
cally localized state at K = 5 and T = 1. For weak in- 
teractions (g = 1) the motion is quasiperiodic with slow 
growth in the number of noncondensed atoms. Strong in- 
teraction (g = 5) destroys the quasiperiodic motion and 
leads to diffusive growth of energy, accompanied with 
exponential growth of noncondensed atoms that clearly 



indicates the instability of the BEC. Notice that the rate 
of growth in energy is much slower than the classical 
diffusion rate, which means that chaos brought back by 
interaction in this quantum system is still much weaker 
than pure classical chaos. 
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Figure 12: Nonlinear effects on dynamically localized states. 
K = 5, T = 1. (a) Plots of average energy E(t) versus the 
number of kicks t, where dash dotted line corresponds to the 
classical diffusion, g = (dash), g = 1 (dot), g = 5 (solid), 
(b) Semilog plot of the mean number of noncondensed atoms 
versus the number of kicks t. g = 1 (dot), g = 5 (solid). 



V. CONCLUSIONS 

We have investigated the complex dynamics of a peri- 
odically kicked Bose-Einstein condensate that is consid- 
ered as a nonlinear generalization of the quantum kicked 
rotor. We demonstrate the transition from the anti- 
resonance to the quantum beating and then to instability 
with increasing many-body interactions, and reveal their 
underlying physical mechanism. The stable quasiperi- 
odic motions for weak interactions, such as anti-resonace 
and quantum beating, have been studied by mapping the 
nonlinear Schrodinger equation to a spin model. The 
transition to instability has been characterized using the 
growth rate of the noncondensed atoms number, which is 
polynomial for stable motion and exponential for chaotic 
motion of the condensate. 

Finally, we emphasize that the results obtained in the 
paper are not limited to BEC and can be directly applied 
to other systems whose dynamics are governed by the 
nonlinear Schrodinger equation. 
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